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A method is presented for obtaining the e-expansion 
for on-shell massless scalar double-box diagram. 

1. There are several reasons why calculation of two-loop diagrams is important. The study 
of two-loop diagrams is required by increasing precision of experiments. Examination of multi- 
jet processes requires analytical calculations of two-loop diagrams. At last, consideration of 
such diagrams helps us to develop and apply new methods of loop calculations. 

One-loop on-shell four-point diagram in cp 4 -theory has been considered in [1]. As 
compared with the one-loop case, the evaluation of two-loop diagrams is technically much 
more complicated and less investigated. Off-shell double-box problem in 4-dimensional space- 
time was solved in [2]. On-shell double-box problem is more complicated than the off-shell one 
because the appearence the infrared divergences. Our consideration will be in the framework of 
dimensional regularization [3]. On-shell double-box diagram has very high degree of 
divergences. Let n = 4+2e be the space-time dimension. Then on-shell double-box diagram 

D(s,t,kf) diverges as £~ 4 . It is evident that the expansion up to the forth degree is a very 
difficult problem. On-shell problem can be solved for example by using various reduction 
methods (the integration-by-part method [4], uniqueness method [5]). Such approach was used 
in [6], but result was too combersome. For this reason it remained the problem to look for 
such scheme, that could lead to more visible result. Such method I'll consider in this paper. 
The talk is about so named Mellin expansion method. At first this method was used for 
multiloop calculations in [7]. With the help of this method I managed to find the visible 
solution for the on-shell double-box case. 

The remainder of the paper is organized as follows. In Section 2 we list notations, some 
most useful formulas and the main idea of calculation. In Section 3 we'll carry out some spade- 
work before to use the Mellin expansion method. In Section 4 using the Mellin expansion 
method we'll obtain the spectral representation for four-point function. In Section 5 we discuss 
the results. 

2. In this section we list notations, some most useful formulas and the main idea of 
calculation. 



D(s,t) = 




Fig.l 

We use the standard notations s = (k l + k 2 ) 2 , t = (k 2 + k 3 ) 2 , k 2 =0. 
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Not to complicate formulae, we omit henceforth the factors that are powers of 2, %, i. These 
factors can be easily restored in the final results. Each line of diagrams carries a power- like 
factor that is pictured as a line with a label a, for example, 



a 
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— in momentum space , 



(x 2 + iO) 



— in coordinate space . 



We'll use constantly the formula for propagator Fourier-transformation from momentum to 
coordinate space: 



1 r(u-a) 1 



(k 2 ) a T(a) {x 2 f~ a 



The main idea of calculations consists in the possibility to reduce in the case k\ = k 2 = 
double-box diagram D(s,t) to three-point diagram (Fig. 2) T(q 2 ,q 2 ,q 2 ): 
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where 

q 2 = (k l + k 2 x) 2 , q\ = (k 3 + k A y) 2 , q\ = txy , x = l-x, y = l-y, 

with the help of equality 
i i 



D(s, t) = j dx j dy T(q 2 , q 2 , q 2 ) 
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(1) 



The three-point diagram is more simple object for calculation and for expansion in 8 than four- 
point diagram and if we will know the e-expansion for T, we'll obtain simply 8-expansion for D 
from this equality. 

3. In this section we'll carry out some preparatory work before using the Mellin-expansion 
method. The reason for this will be clear later. Using the uniqueness method we present the 
three-point diagram (here it is written in coordinate space) in the following form; 
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We begin our calculation from the usual Feynman a-representation for diagram of Fig. 3 in 
coordinate representation. Then we transform this a-representation to spectral representation 
for Fig. 4. 
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Fig. 3 Fig. 4 



4. In this section we obtain at first spectral representation for diagram of Fig. 3, then using 
the Mellin expansion method we'll reduce it to spectral representation for four-point function. 
Spectral representation for diagram of Fig. 3 has the form 
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This expression can be presented in the form of two-fold Mellin integral [8]. After 
multiplying this Mellin expansion by the propagators that correspond to formula (3) and after 
transforming the result to momentum representation we obtain for two-loop three-point 
diagram of Fig.4 in momentum representation the following a -representation 



fi-rj 2/1 \ afafod 2e (p(oc ( ) 



where 



\ x J {a l xz + a 2 z + a 3 } 



l+e 



4 



This representation is very similar to ordinary a-representation for one-loop three-point 
function, but the representation for two-loop case include weght- function cp(a ■ ) . 

After transformations that take into account equality (1) we obtain for four-point function: 

1 i oo oo *'ee2e °° 

D(s,t) => jdxjdyjJZjdZ V; J^x x-Hl + x) 1 " 

{(7 2 +^ + q {^ + «?1^ + «?2%| o 
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x ■ ■ 



: f<fe 

C {(l + xV+S + z} 
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At this stage we can see why the previous transformations were done. For the form (3) we can 
carry out integration over y. Then 

and taking into account, that q 2 = sx we obtain (here L = t I s) 
D(s,t) =>—rwrJ, 



3-2E 
S 



J i {1 + T + ^}{1+T + ?} 
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{xC + (l+T)(x + ^ + Q} o {XZ + X + Z + ^} 
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This spectral representation is a convenient starting point for both analytical and numerical 
calculations. For instance, in a separate publication I will obtain one of the irreducible diagrams 
that appeared in [10], which corresponds to L = -1 as well as a much more cumbersome 
analytical result for arbitrary L as a quadratic form of Li n -functions. 



Conclusions. In this paper we have considered an approach to the calculation of on- shell 
two-loop scalar double-box diagram. This approach is based on the following tools: (i) the 
Feynman parametric representation, (ii) the uniqueness method, (iii) Mellin-Barnes contour 

integrals. We considered only scalar diagrams (corresponding to massless cp 3 -theory) because 
expressions occuring in realistic calculations can be reduced to such scalar integrals. We have 
obtained the analytical formula (6), that can be used for calculations in various concrete cases. 

I am grateful to F.V.Tkachov for his interest in this work, useful discussions and help. The 
reseach was supported in part by the Russian Fund for Basic Reseach (grant 95-02-05794). 
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